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Abstract 

A simplified Heisenberg spin model is studied in order to examine 
the idea of decoherence in closed quantum systems. For this purpose, 
we present a quantifiable definition to quantum coherence S, and dis- 
cuss in some detail a general coherence theory and its elementary re- 
sults. As expected, decoherence is understood as a statistical process 
that is caused by the dynamics of the system, similar to the growth of 
entropy. It appears that coherence is an important measure that helps 
to understand quantum properties of a system, e.g., the decoherence 
time can be derived from the coherence function S(i), but not from 
the entropy dynamics. Moreover, the concept of decoherence time is 
applicable in closed and finite systems. However, in most cases, the 
decay of off-diagonal elements differs from the usual exp(— t/r^) be- 
haviour. For concreteness, we report the form of decoherence time Td 
in a finite Heisenberg model with respect to the number of particles 
N, density n p , spatial dimension D and e in a rj/r e -type of potential. 

1 Introduction 

Decoherence is widely accepted as an explanation of how quantum corre- 
lations are damped out to make physical systems effectively classical, and 
thus it provides a satisfactory answer to the dilemma of Schrodinger's poor 
cat [l]. Open quantum systems with an infinite environment have been 
studied in great detail in various environment-induced decoherence schemes 
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El HI El El [7j using master equations of reduced density matrices. How- 
ever, the "cosmological" aspect of physical reality motivates studies in closed 
quantum systems since, strictly speaking, the universe has no environment 
[31 SI El [9] and yet decoherence is observed in our universe (TQl [TT] . The most 
interesting questions are: How coherence in the universe should be modelled 
and understood? How reasonable it is to talk about decoherence in closed 
systems? There are three main lines in studying closed quantum systems: 
decoherent histories approach [12l EH HH E3 HE], self-induced decoherence 
scheme pH HE! EH EHJ EE] and density matrix formalism. Density matrix 
formalism has also been used to study open and finite quantum systems 
[El EH EEl [23]. 

Decoherent (or consistent) histories approach uses the frame of the many 
histories interpretation of quantum mechanics. The source of decoherence is 
projections along histories that produce effective coarse-graining, and thus 
decoherence. Intuitively, histories are analogous to paths in path integral for- 
malism, and meaningful histories should obey consistency conditions. The 
criticism against decoherent histories approach mainly consists of the fol- 
lowing arguments: (1) The notion of "decoherence" in decoherent histories 
approach is not equivalent to "decoherence" in density matrix formalism [13J. 
(2) Philosophically, decoherent histories approach seem to violate one of the 
main principles of quantum physics, i.e., in principle, the wave function is 
the complete description of the state of the system. Histories seem to be the 
complete description and the wave function incomplete |24| . (3) In decoher- 
ent histories approach, decoherence effects are dependent on the chosen set 
of histories [251 [26] • The choice of the set of histories is arbitrary and results 
in contrary propositions [27], i.e., physical phenomena that are predicted (or 
retrodicted) using the consistent histories approach depend on the chosen set 
of histories. In Ref. [28] a peculiar feature of decoherent histories is shown: 
in a universe in which decoherent histories approach is valid, the objective 
properties are not objective at all, but dependent on chosen set (or family) 
of histories. This is a version of the preferred basis problem [8l 129]. 

In the self-induced decoherence scheme, decoherence is created within the 
undivided system. There is no need for environment as the source of decoher- 
ence, unlike in the environment-induced decoherence (in open systems). De- 
coherence appears because quantum correlations of expectation value (M) in 
continuous energy spectrum are proportional to / dE f dE'M(E, E')e~ l ^ E ~ E 
and, according to Riemann-Lebesque theorem [30], it vanishes when t — > oo 
if the integrand is L 1 integrable, and thus results in a diagonalisation of 
expectation value. There are three crucial assumptions: (1) continuous en- 
ergy spectrum (i.e., the system is infinite), (2) the coefficients in the energy 
eigenbasis are integrable functions of energy (for finite systems this condition 
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naturally fails), and (3) the time scale is infinite [9]. As analysed in Ref. [9], 
for finite systems in long (but finite) time the diagonalisation does not occur. 
However, the critique against self-induced decoherence scheme presented in 
Ref. [9j is based on open and finite system analysis that is fundamentally 
different from closed and infinite system, as demonstrated in Section 12 . 1L 
Still, the definition of coherence in self-induced decoherence scheme lacks 
generality since it is applicable only to infinite systems. 

The density matrix formalism is a natural way to study quantum systems, 
so it should be useful also in decoherence studies in closed systems. Moreover, 
it does not have the same drawbacks as the two other main lines to study 
closed systems. However, the concept of coherence is not clear at all (neither 
in density matrix formalism nor in any other formalism). A well established 
coherence theory using density matrices should help clarify the concept of 
decoherence in closed systems. Coherence theory and conceptual problems 
of decoherence (especially in closed and finite systems) are interesting and 
wide research topics alone, but related to this study, the following elementary 
results are important. 

1. There is a basis- vector-set-independent way to quantify coherence uniquely. 

2. The possibility of recoherence does not mean that there is no decoher- 
ence. Decoherence is a dynamical phenomenon that reduces coherence, 
and hence, recoherence is a dynamical phenomenon that increases co- 
herence. All finite quantum systems may experience recoherence. 

3. There are two different coherence types (similar to the different entropy 
types): idealistic and realistic coherence. 

• The idealistic coherence scheme can be applied only by those ob- 
servers who do not interact with the universe, and who know the 
wave function of the universe and its time evolution. In a closed 
system, there is no idealistic decoherence, since the wave function 
of the universe always remains pure. 

• The realistic coherence scheme is the internal view of the universe 
calculated from the wave function of the universe. It describes 
the events as the real observers that are totally correlated with 
the universe perceive them. To acquire this realistic viewpoint, 
one should make an effective theory from the wave function of 
the universe, e.g., to use reduced density matrices. This is often 
referred to as coarse-graining. 
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A short analysis of the results is presented in Section El In this study, 
the word "decoherence" refers to the concept of realistic decoherence. 
Exceptions are mentioned. 

Idealistic decoherence is often referred to with the word "decoherence", since, 
in open and infinite system studies both types of decoherence behave sim- 
ilarly. The principal difference between idealistic and realistic decoherence 
can be seen only in closed systems. 

The first main purpose of our study is to define coherence uniquely and 
to analyse carefully the consequences of coherence defined as 



where A max is the maximum eigenvalue of the density matrix p, and N is the 
dimension of the Hilbert space of the corresponding quantum system. The 
definition satisfies the requirements for coherence and is useful in, e.g., calcu- 
lating decoherence times of quantum systems. However, it is not an extensive 
measure like entropy, but intensive and probability-like. Both entropy and 
coherence are important concepts in understanding different properties of a 
quantum system. 

Second, we apply the above presented coherence function to study a 
Heisenberg spin model in order to clarify the idea of decoherence in closed 
systems. We focus on the single particle coherence, and thereby calculate 
the time evolution of the (single particle) coherence of the system. This 
coarse-graining is similar to the one used in Gibbs coarse-grained entropy 
calculations [HI [32]. It is demonstrated that S is a valid and useful measure 
for coherence and coherence dynamics of the simple closed and finite system 
under the study. Moreover, we determine how the decoherence time depends 
on the system variables (particle number N, density n p , potential and dimen- 
sion D of the system). The behaviour of decoherence is calculated for i]/r e 
potentials in various dimensions D. The results are interesting: the usual 
coherence decay profile exp(— t/r d ) occurs only in special cases. Generally 
the decay follows 
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where f(D, e) = 1.97 (1 - 0.93 exp (-0.65D L35 e- L68 )), and 
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is the decoherence time. 
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This paper consists of six main sections. First, we review a relevant part 
of coherence theory in Section El The spin model is introduced in Section 
El In Section [4] we present theoretically the dynamics of a simple initial 
state, and apply the coherence function E(t) as the measure of the realistic 
coherence of the whole system. The structure of numerical simulations is 
presented in Section El and Section [6] explains the results of the numerical 
simulations along with an analytical reality-check of basis independence of 
our model. Section [7] is for discussion. 

2 Introduction to general coherence theory 

In a scientific dictionary, coherence is defined as "correlation between the 
phases of two or more waves, so that interference effects may be produced 
between them, or a correlation between the phases of part of a single wave" 
and "property of moving in unison" [33]. The amount of coherence of waves 
depends on how correlated they are; how close the wave lengths, phases and 
amplitudes are to each other. Thus, coherence can generally be defined as 
"correlation of wave-like entities" and it enables interference. Quantum phys- 
ical objects are generally described using wave functions or state vectors, 
which are "wave-like entities". Thus, quantum coherence could be under- 
stood as "a measure of the strength of quantum correlations". Within the 
framework of classical correlations, quantum correlation phenomena cannot 
be explained. Evidently, entanglement and superposition states in general 
can produce observable non-classical correlations [34]. In the density matrix 
formalism, the fingerprints of quantum correlations are the off-diagonal ele- 
ments of the density matrix. From now on, in our article the word "coherence" 
mostly refers to quantum coherence. 

Let p be a density matrix of a closed system in the basis set A = {\cti)}, % = 
1 . . . N. The basis vectors span a ^-dimensional Hilbert space H. The equiv- 
alent description about the (closed) system can be given with an arbitrary 
basis set B = if the basis set spans the same Hilbert space H. The 

quantum dynamics of the system does not depend on the chosen basis vec- 
tors. 

Let us now consider quantum correlations (coherence) in a basis set, say 
A. The word "quantum coherence" refers to the off-diagonal elements of the 
density matrix p, but the height and location of off-diagonal elements (and 
diagonal elements) depend on the chosen basis set. Therefore, the "pointer 
basis" has been invented: it is the basis set in which the quantum dynamics 
diagonalises the (reduced) density matrix [HI HJ E2J EH]. If one wants to 
describe the whole closed quantum system without using tricks to make it 
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effectively open, one cannot use the pointer basis since the unitary dynamics 
does not diagonalise the (irreduced) density matrix in any basis. 

However, the question "How much coherence is there in a given density 
matrix of a particular closed system?" is a relevant one, and the answer 
should not depend on the chosen basis set. In a closed system, the quantity 
called "coherence" of the total density matrix p should be a conserved quantity 
because of the previous reasoning. According to the Noether's principle [36j, 
the "conservation law of coherence" is a result of a particular symmetry of 
the system. The best candidate for that symmetry is the invariance of the 
density matrix under the change of the basis vector set. If one wishes to 
define "coherence", the definition should obey these constraints. 

The density matrix of a totally coherent system of N possible states 
described with, say basis set A = {|aj}},i — 1...N, 
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in the eigenvalue basis B = {|A)} that is always diagonal. The transfor- 
mation is obtained by unitary transformation U so that ps = U pjJJ^ x . By 
previous reasoning, these density matrices should contain the same amount 
of coherence. And, because quantum correlations refer to correlations that 
cannot be explained by classical correlations, they should be measurable - 
at least in the idealistic sense. The totally classical density matrix 
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(6) 



is diagonal in any basis. Let us consider following measurement scheme: we 
are free to measure an observable O in any kind of basis. The probability to 
observe "a matter of fact" \0) is 



Prob(p,0) = Tr(pO) = Tr(p\0)(0\). 



(7) 
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While measuring pa, the probability Prob(p^,0) is maximised when \0) = 
-7= (1 ■ ■ -1)^, and thus Prob(pA,0) = jfY^jPij — 1- For the totally clas- 
sical density matrix pc, the probability for the same measurement out- 
come is Prob(pc,0) = jfYnj Pi,j — jf- Thus, maximum quantum corre- 
lations (coherence) are causing a 1 - | increase in the probability of the 
measurement outcome. Naturally, the same probabilities of measurement 
outcomes happen in the eigenvalue basis, so that the measured state is 
\0) = U\0) = (1 ■ • • 0) B . The eigenvalue basis is useful, since the max- 
imum eigenvalue A max is the probability for a state (of the measurement 
outcome) to be in its most probable state, i.e., to demonstrate maximally 
the quantum nature of a particular system. 

Thus, we define coherence 5 of the density matrix p by 

= N-i (" )vmax " 3 ~ jy) ' ( 8 ) 

where N is the dimension of the density matrix p. Subtraction of the term 
1/N, which represents the classical probability, is a background correction 
that sets the minimum coherence to zero. The prefactor normalises the 
maximum coherence to 1, since a maximally coherent system would yield a 
probability-optimised measurement outcome with probability of 1. The value 
of S describes how strongly the quantum coherence influences the behaviour 
of a particular system as whole. Since the value of coherence S depends only 
on the dimension of the Hilbert space (N) and the maximum eigenvalue A max 
of the density matrix p, it clearly is a basis-set-independent measure. 

2.1 About coherence dynamics 

There are some important remarks concerning coherence dynamics of closed 
quantum systems. At first, let us consider a classical example of, say 10 9 , 
longcase pendulum clocks that do not interact with each other. They are 
similar to each other except that their periods are random, due to random 
effective lengths of pendulums. At the instant of time t the phase of all pen- 
dulums is the same, i.e., the pendulums are at the lowest possible point and 
swinging leftwards. Classically thinking, this initial setting is very coherent 
in the phases of pendulums. After a short time, at the moment t\ the phase 
coherence is totally lost. However, at the moment tp the phase coherence has 
revived - pendulums are at the lowest possible point and swinging leftwards 
- since the number of pendulums is finite and they all have fixed periods. 
The revival time tp is known as Poincare recurrence time, and is the longest 
possible timescale in the system. In an infinite time, a finite system recurs 
arbitrarily close to its initial state (or any state) arbitrarily many times. [38] 
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Let us consider now a Gaussian wave packet that is trapped in an anhar- 
monic potential U(r) = \muj 2 {r — r ) 2 + x(r — r ) 3 that has small anhar- 
monicity x- At first, the wave packet oscillates coherently in the potential. 
After an intermediate time the wave packet has lost its form and shape, i.e., 
has lost its coherence. However, after a very long time the wave packet re- 
forms again to Gaussian shape if the system consists of a finite amount of 
bound states. For more details, see Ref. [39] and references therein. The be- 
haviour of coherence can be studied via the time evolution of the expectation 
value of position (x(t)). The strong oscillations of the expectation value of 
the position present in the beginning, dampen as the system evolves in time, 
resulting in a fluctuating mean. 

The third interesting example is interacting bosons in a boson field so 
that the group of bosons forms the system and the boson field forms the en- 
vironment. The standard procedure to study the problem is to trace out the 
boson field from the complete description density matrix in order to acquire 
the master equation for the density matrix of the system. The system is an 
open one. See, e.g., Refs. [It)]: Ch. 14, [HJ: Ch. 6.1, [12]: Ch. 1, [43]: Ch. 
3.3. At first the boson field is treated to have a finite (but large) number 
of degrees of freedom, but when correlation functions should be calculated, 
the number of degrees of freedom of the boson field is assumed to be infinite 
(justified by the fact that it is very large compared to the number of degrees 
of freedom of the system), so that the sums can be converted to calculable 
integrals. These kinds of environments are accepted sources of decoherence, 
resulting in a monotonic decay of coherence that has characteristic decoher- 
ence time [7]. 

In all these examples the concept "coherence" (classical or quantum) is 
reasonably applicable, but the concept "decoherence" is commonly associ- 
ated only with the third example because that example describes a net loss 
of coherence out of the system. The fundamental difference between the sys- 
tems of the first two examples and the last is that the system of the last 
example is an open system and an infinite system. The net loss of coherence 
requires both "openness" and "infiniteness" - openness because an environ- 
ment is needed for coherence dumping, and infiniteness because from a finite 
environment all coherence will eventually come back to the system. Thus, 
the observed net loss of coherence in some open quantum system studies is 
the result of an approximation that assumes infinity because of the need of 
calculating integrals instead of sums. It is an approximation that changes the 
fundamental behaviour of the system: it results in a monotonic decay of co- 
herence with no revivals. Usually, in the fundamental level, approximations 
neither increase accuracy nor enlarge the area of application. Moreover, it is 
not even clear that our universe is infinite. The Bekenstein bound [44] sets 
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the upper limit to information (related to the degrees of freedom) in finite 
volume, and estimated (quantum) computational capacity of our universe is 
finite [45]. 

The longcase pendulum clock example demonstrates the difference be- 
tween a finite closed system and an infinite closed system: if one increases 
the amount of clocks, Poincare recurrence time either increases or remains 
the same. There is no increase of recurrence time if the period of the added 
clock multiplied by an integer is the period of the existing clock system. If 
periods are random real numbers (it does not matter if the interval is fixed), 
for infinite amount of clocks the Poincare recurrence time is infinite, i.e., there 
is no recurrence. Even if the system is closed, the decay of phase coherence is 
monotonic. If one considers the second example system with infinite setting, 
this feature is even more striking. The system is closed, so that there is no net 
loss of coherence out of the system. Infiniteness means there is no Poincare 
recurrences, because the initially very narrowly concentrated finite amount of 
coherence spreads evenly to infinitely small amounts of correlations between 
all the infinite states - and it never finds its way back. 

While closed and finite systems differ fundamentally from open and infi- 
nite systems, it appears that the term "decoherence" is an accurate descrip- 
tion for what happened in these thought experiments. The word "decoher- 
ence" stands for a dynamical phenomenon that reduces coherence. The word 
"decoherence" does not even suggest that this reduction of coherence should 
be permanent. Finite systems may experience both decoherence and reco- 
herence. Recoherence is a dynamical phenomenon that increases coherence. 
In density matrix formalism, decoherence is commonly understood as the de- 
cay of off-diagonal elements of the density matrix. More precisely, in general 
case, decoherence of a particular density matrix p is the decay of the function 
S(p). 

2.2 Realistic and idealistic standpoints 

Similar kind of questions about the behaviour of closed systems have been 
a leading topic of debate over the entropy and universe more than hundred 
years. Here we will give only a brief introduction to this topic. A more 
interested reader is encouraged to study the vast literature covering the issue, 
starting with [3H E21 HQ HZ] HSl \M ED] and references therein. 

Entropy is an "extensive measure depending on the number of possi- 
ble (micro)states and corresponding probabilities". The probabilities of mi- 
crostates are used as weights in Gibbs entropy Sq = — Y^fLi Pi hipj. Gibbsian 
entropy is often titled as "statistical entropy", as well as its quantum equiv- 
alent, the von Neumann entropy Sn = — Trplnp. Entropy obeys the same 
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"independence of the basis vector set" principle as coherence. 

The second "law" of thermodynamics states that the entropy of a closed 
system cannot decrease. However, the underlying laws of physics are re- 
versible in time, and if the closed system is finite, then it should experience 
Poincare recurrences, which means a decrease of entropy at some moments. 
Another entropy problem that has relevance to coherence studies was the 
observation that the entropy of a closed system (universe) studied as whole 
was a constant of motion while inside the universe it seems that entropy is 
increasing basically everywhere. 

The irreversible entropy dynamics that is associated with the second "law" 
was a result of an approximation in deriving the Boltzmannian H-theorem 
[Tfl I48|. However, even Boltzmann himself noticed that irreversible entropy 
dynamics is not a lawlike theory but merely a statistical theory [51] that 
allows the decrease of entropy. The second entropy problem was partially 
solved by introducing a method to coarse-grain the complete system de- 
scription so that the statistical behaviour of entropy remained unchanged 
in closed and finite systems. In Gibbsian coarse-graining, probabilities of 
infinitely small phase-space elements are averaged in a finite phase-space 
element, and then entropies of all averaged probabilities are summed to ob- 
tain the entropy of a closed system [3U [32]. However, this did not end the 
philosophical discussion in which the dissatisfaction concentrated basically 
in two kind of arguments: "what is the 'right' coarse-graining to choose from 
the many possibilities (i.e., how to divide the system into subsystems)" and 
"while the 'true' entropy of closed system is a constant, what does it mean 
that coarse-grained entropy changes in time". 

Let us now introduce two notations for entropy and coherence of a quan- 
tum system described with the density matrix p: idealistic entropy = 
— Trplnp, idealistic coherence S ic j = H(p) and realistic entropy SVe = — X^Trpdnp;, 
realistic coherence E re = iV -1 £^ S(pO , where pi is the density matrix of a 
subsystem. For any closed system, idealistic entropy and coherence are con- 
stants of motion. In calculating idealistic entropy and coherence, everything 
inside the closed system is treated quantum mechanically, and if the universe 
is under consideration, it still demands that measurements, measurement 
apparatuses and even observers are treated quantum mechanically. Mea- 
surements are expressed as interactions between measurement apparatuses 
and measured systems and so on. Clearly this is possible only for an idealistic 
observer, say, a Clever Chinchilla. She does not interact with the universe 
(closed system), and after having guessed the correct wave function and its 
time evolution, studies the universe only as an academic example in her brain 
(or with a pen and paper) without ever interacting with it. 

The realistic standpoint describes what the real observers inside the uni- 
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verse observe. These observers are totally correlated with the universe and 
can only in their wildest imagination reach the Clever Chinchilla's view of 
the universe - as we do in this article, when we refer to idealistic quanti- 
ties. If the Clever Chinchilla wishes to study the universe of her academic 
example in order to figure out what the universe looks like when viewed from 
inside, she coarse-grains her perfect model and obtains realistic quantities. 
Naturally, idealistic and realistic quantities behave differently, but there is no 
contradiction between them - if the coarse-graining is done right (note that 
the Clever Chinchilla is so clever that she never makes mistakes), she ob- 
tains as accurate a description as possible of the different possibilities of the 
universe as the Real Observer would see it. Because quantum processes are 
truly random, the Clever Chinchilla only acquires all possible measurement 
outcomes weighted with their probabilities. 

Both realistic and idealistic quantities describe the same system, but from 
different standpoints. Entropies are related via the quantity I = S TC — Sid, 
which is called mutual information [52]. Mutual information expresses the 
amount of correlations between the (arbitrarily chosen) subsystems. Close 
relations between entropy and coherence suggest that there exists also a 
quantity E = H ic i — H re that quantifies the degree of mutual entanglement, 
i.e., how entangled the subsystems are with each other. This quantification 
scheme for entanglement is applicable for all quantum systems in density 
matrix formalism. 

The perfect description of the system is impossible inside the closed sys- 
tem, even in principle. The fundamental logical difference between idealistic 
and realistic standpoints is similar to Godel's incompleteness theorem [54], 
which states roughly that in a richer logical system than the first order predi- 
cate logics (the universe) that has propositions concerning its own consistency 
(Real Observer, with the help of measurements, searches the perfect descrip- 
tion of the universe) there exist propositions <fi of the form of cannot be 
proved true". To prove <fi requires a higher order metatheory (Clever Chin- 
chilla). Thus, the closed system cannot be completely described within the 
system, and the complete and incomplete descriptions are different. 

Decoherence is an observed phenomenon [10, CE] in the universe. Thus, 
the decay of coherence should refer to the quantity S re - even if the whole 
universe obeys unitary dynamics. There is, however, one problem without 
a proper answer: how to make coarse-graining correctly - how to divide a 
closed system into subsystems. In this study, we coarse-grain to the smallest 
level of individual particles. A more detailed analysis of coherence theory and 
problems related to defining decoherence will be published elsewhere [55] . 



11 



2.3 Entropy and coherence 

Both entropy and coherence are important concepts that help us to classify 
and understand properties of quantum systems. Entropy is an extensive 
measure, while coherence is intensive and behaves like probability [56]. Let us 
illustrate some properties of entropy compared with properties of coherence 
by considering two identical uncorrelated spin-| systems that both have, say, 



Pi = ( ? f ]■ The eigenvalues are A± = | ± |. Thus, both systems have 



coherence Si = | and entropy Si ~ 0.56. The system p i+2 = p\ ® p 2 has 

eigenvalues A = ^A^A^A^, and thus the entropy is S\ + 2 ~ 1.12 and 
entropy per particle is S norm ~ 0.56, and coherence is S 1+2 = ^. Moreover, 
constructing a greater system from uncorrelated smaller systems causes a 
dramatical drop in coherence, since the coherence of N joint uncorrelated 
systems is proportional to Hf A max ,/. 

However, the above mentioned "anomaly" in coherence is not surprising at 
all, if one remembers that coherence is closely related to probability, that is, 
to the maximum possible probability of the system. If one has the probability 
of 3/4 to measure the maximum probability state of a two-level system, the 
probability of maximum probability state for two such systems is (3/4) 2 . 

If one considers unitary dynamics, both (idealistic) coherence and (ideal- 
istic) entropy of closed system are constants of motion that give constraints 
to the system. In the system of pi +2 the transition to the system that 
has eigenvalues ^ A ^ would be "legal" from coherence-point-of-view, but 
it would change the entropy S1+2 — 1.12 — > 0.69, and while eigenvalues 
A ~ 0.5 A 0.32333 A 0.125 A 0.05167 are possible from entropy-point-of-view, 
they would result in change of coherence of the system. Coherence has signif- 
icant importance if one wishes to compare quantum properties of two systems 
that, e.g., have the same entropy. Of those, the system that has the greater 
coherence also has more quantum correlations, and would be more useful to, 
e.g., quantum computation. Entropy does not tell us the whole story. 



We have decided to study the Heisenberg spin model because it is simple 
enough to solve, and yet complicated enough to simulate the properties of 
real quantum systems. Coupled spin systems are interesting from a quantum 
computational point of view, too. Our system, N interacting particles fixed 
in a space, has no environment and, in that sense, the system forms a closed 
quantum universe. The particles are spin-| particles, and the interaction be- 




3 The model 
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tween them is due to their spin-z component (analogous to the Ising model). 
When there is no coupling with the environment (i.e., no outer environment), 
the two spin states have the same energy, which is taken to be zero. Zurek 
[22], Omnes [23] and Schlosshauer [HI [9] have considered a similar, but sim- 
pler model in order to study the decay of off-diagonal elements (i.e., quantum 
correlations) of a reduced density matrix. They labeled one particle as the 
system, and the others as the environment [59], but, in this case we study 
the particle system whole. 
The interaction Hamiltonian, 

N-l N N 

H = hY, E 9a< ® < II ( 9 ) 

j=l i=j+l k=l,^j,i 

describes the dynamics of the system. The interaction matrix G, where 
9ij = 9ji> gi yes the interaction strength between particles i and j. The 
interaction strength arises from the potential V; but, for formal calculations 
there is no need to know more about it, because particles are doomed to stay 
in one place. Fixing the positions of the particles is a justified assumption 
in decoherence studies since, in most cases, the decoherence time scale is 
the shortest time scale [3], at least shorter than the time scale of particle 
motion. In numerical simulations, only potentials of the type V = r)/r e are 
considered. 



4 Theoretical calculations 



Let us now consider a simple case in order to present our method. The task 
is to solve the time evolution of the realistic coherence S re (t) of a closed and 
finite system that is initially in a product state of superposition states 



N 



i*(o)> = n +h\-k)) 



k=l 



where a's and b's are normalised probability amplitudes |afc| 2 + \bk 
all k = 1, . . . , N. The Schrodinger equation, 

mMt)) = H\*(t)), 



(10) 
| 2 = 1 for 

(11) 



gives the dynamics of the system, and with the given initial condition of 
equation ( TTOl) one gets the time dependence 



|*(t)) = exp 



N-l N 

j=i i=j+i 



A' 



II ® Kl+fc) + h\-k)) 

k=l 



(12) 
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The fate of the I th particle is solved by tracing over other particles, i.e., 
degrees of freedom, 

Pi = Tri,...,^ p , (13) 

where p = \^(t)){9(t)\. This is a crucial step. With Gibbsian coarse-graining 
[BTl [32] we make an effective theory of our particle system by tracing over the 
"uninteresting" particles (that form an effective environment to the particular 
particle of our interest), as in the mean field approximation. The net effect 
of traced-out particles is described in a simpler form and with less degrees of 
freedom. 

We thus have 



Pi = \ai\ 2 \+i)(+i\ + \bi\ 2 \-i){-i 
+ 



N 

aft II (k| 2 e~ i2ff!fc * + |^| 2 e i2to ') |+,)(-,| + h. c. 



k=l,^l 



(14) 



It is interesting that the result of Eq. (1141) is the same as in Ref. [22], if one 
drops off the index I. In Ref. |22j only an interactionless environment has 
been considered, but our model counts all the interactions between particles. 
Let us make our notation a bit lighter by denoting 

JV 

Zl = ai b1 I] (k|V i2 ^ + |& fe |V 2 ^). (15) 

k=l,^l 

This z\ (or its complex conjugate z*) describes the fate of the off-diagonal 
elements of I th particle. The eigenvalues of the reduced density matrix p\ are 

1 1 



A w = -±- > /l-4(|o,|W-|^)| a ). (16) 

Now the time evolution of single particle coherences S/ sing i e = S(/?j) and 
the realistic coherence S re (t) = A^ -1 Ya=\ ^-z.singie of the closed system can be 
evaluated. Inserting maximum eigenvalues of Eq. f TT6l l into the definition of 
Eq. (jHJ) results in 

„ . , 1 M / 1 \ 1 "2 (\ 1 r — — - . ; 1\ 

=" (( hv g M=i ( x -' - m) = n g 1 (2 + 2iA-«<W , W-W>l ! '>-5) 
1 N I 

=T7 E V / l-4(|az| 2 |6z| 2 - kz(t)| 2 ), (17) 
iV i=i 

where is the number of particles and M is the dimension of the reduced 
density matrix. 
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5 Numerical simulations 



The main aim of our numerical simulations is to demonstrate the usefulness 
of the definition of coherence presented in Section El The secondary aim is 
to study coherence behaviour of a simple closed and finite system, includ- 
ing Poincare recurrences. The third aim is to calculate the dependence of 
decoherence time of relevant system variables, i.e., = r^(N, n p , rj, e, D), 
where iV is the number of particles, n p the particle density, 77 the strength 
of interaction of rj/r e -potential, and D the spatial dimension of the system. 
The usual coherence decay in infinite systems has the form of exp(— t/r d ), 
which is used to define the decoherence time [31 H]. It appears that the decay 
profile of the form of exp ( — (t/T^)^ '^) is more accurate than the usual one. 



The starting point of simulations is a D-dimensional "box" whose volume 
is l D . N particles are placed randomly in this box. These particles are in 
fixed places, and they interact with each other according to the Hamiltonian 
([2]). The size of the box is related to the number of particles N and the 
particle density n p by 

fN\ 1/D 

' = UJ ' (18) 

Without lack of generality and to ease numerics, we choose to study an initial 
state that contains complete superpositions, so = bk = l/V%- Thus, 
the closed and finite system has idealistic coherence = 1 and idealistic 
entropy = 0. The interactions between particles have the form of gik = 
T]/rf k , where is the distance between / th and k th particle. Inserting these 
conditions to Eq. (fTTj) results in the realistic coherence 

s»(*)=4e n cos (?V ( i9 ) 

iV 1=1 k=l,& V T lk ) 

Second, we use a least squares fit on S re (t) using a function 

£(t) = (1 - c)e- (t/ * d)C + c (20) 

that allows fluctuations around the average level c of S re (t). The average 
level c is the expectation value of H re (t) in the interval [t\, £2], where <^ t\ 
and r d <t 2 -ti: 

c = (5 re ) = f 2 E Ie (t)dt. (21) 

C2 — 1\ Jt\ 



The sum of squares error \ f° r the function of Eq. (1201 describes the ac- 
curacy of the fit. On the last phase of numerical analysis, it will be used 
to obtain relative weights for data points. The data of one run consists of 
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fitting parameters t d (represents the estimate of decoherence time), C (expo- 
nent that characterises the coherence decay of r~ e -dependent potential in a 
certain dimension D), c (the fluctuation level of coherence) and x- Poincare 
recurrence times Tp are calculated directly from gikS. 

The procedure is repeated U times. For U = 100 we get a good statistical 
estimate of the behaviour of the system with particular parameter values 
(Ni,€i, Di). The dependencies about rj and n p are included in the unit of 
"simulation time", i.e., T = r]n e J D t. The parameter space that is studied 
thoroughly is A G [20 . . . 100] A D G [1 ... 4] A e G [1 ... 2] and also some 
additional points and slices from greater A-values, most notably with A = 
200: D E [1 ... 10] A e G [1...10]. 

6 Results 

Numerical simulations clearly show that the proper coherence decay profile 
in our model is = (1 — c)e~^/' d ^ +c. Figured] demonstrates the accuracy 
of the fit. The proposed decay profile fits without significant fluctua- 
tions or deviations to simulation data if D + 1 > ^/A/200exp (e/2). The 
main factor for the observed feature is the behavi our of the distance mea- 
sure of two coordinate points I = — x'i) 2 ', it allows more particles 
within a fixed distance from the reference particle as the amount of spatial 
dimensions D increases. Another factor is that as e increases, the range of 
potential shortens and thus less particles within a fixed distance contribute 
to coherence decay. Figure [2] presents the behaviour of expectation value 
of fluctuations and average of weights in respect to dimension D and the 
potential term e. The minimum fluctuation level behaves as ~ lO^^ 5 . Fluc- 
tuations begin to give a strong contribution if D + 1 < ^ A/200 exp (e/2), so 
that fluctuation level with D = 1 and e = 10 is ~ 0.3 in which case fluctua- 
tions of coherence are quite dominant and thus that kind of system possesses 
considerable amount of its quantum behaviour despite the decoherence. As 
stated above, it is not surprising at all, since the short range of the poten- 
tial isolates particles practically from all other particles except their nearest 
neighbours. However, the ~ v^A behaviour for the low-level fluctuation limit 
results in the interesting fact that even systems with e = 1 -type of potentials 
do not lay on lowest coherence fluctuation level if the particle number is large 
enough. It seems that this result allows considerable coherence fluctuations 
for macroscopic objects, but the following facts should also be taken into 
account before claiming that these fluctuations will play significant role: (1) 
The floor of the fluctuation level drops as ~ 10 _7V//5 , so that the net effect 
of coherence fluctuations may still get smaller even if the system parameters 
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do not lay on the diminishing bottom area. (2) The physical effect on the 
particular physical object is a sum of all interactions; including also those 
interactions that are on the bottom fluctuation level. 

The power C of the exponential decay profile seems to have the upper 
limit C < 2. As demonstrated in Figure [U the usual coherence decay occurs 
if D = e = 1. The functional forms of decoherence time and exponent C 
are analysed in detail in Section 16.31 The sum of squared errors y 2 is used 
as the inverse of weight Wi = Xi 2 ■ 

6.1 About interactions and decoherence in different ba- 
sis sets 

Let us study decoherence of the presented well-defined closed and finite model 
in a different basis set. So far, the study has considered the problem only 
in the a^-basis, but, the perfect superposition state in <r 2 -basis is a well- 
defined eigenstate in <r x -basis (spin up). According to the physical intuition, 
the results of nature should not depend on which way they are theoretically 
described. Thus, we should do a reality-check and study what happens to 
the density matrix in a general spin basis. 

One can obtain any spin-basis via unitary transformation 

cos 6* e~ l(f) sm9 \ , . 

e** sin 9 -cos 9 J ' ^ > 

The initial state of Eq. (fTOl) . with = b k = 4, is in the general spin-basis 

N 

l*(0)W= II® (( cosd + e- i<f, sm6)\+ k ) e>4 ,+ (-cos0 + e^sin0)|- fe }^) . 

k=l 

(23) 

If one chooses 9 = 7r/4 and = 0, one gets cr^-basis which is a product 
state of eigenstates with the probability of 1. This initial step of analysis 
illustrates some problems of the decoherent histories approach: how would a 
pure state with the probability of 1 experience decoherence? One can always 
find a basis set in which a superposition state appearing in another basis set 
is an eigenstate. 

One should remember that decoherence is a dynamical process, and thus 
the interaction Hamiltonian is an equally important part of the problem as 
the initial state. In the present case, the Hamiltonian of Eq. is in ag^- 
basis 

N-l N N 

H = hJ2 E 9ij<4,* ® a k<P II ®!*> (24) 

j=l i=j+l k=l,jt=j,i 
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where 

( cos 2 9 -sin 2 6 e - ** sin 20 \ , . 

^ ~ sin 20 sin 2 - cos 2 " (25j 



Thus, the reduced density matrix in a general spin-basis is 



| + sin cos (e -i *2i (t) + e^zf (t )) e" i2 ^ sin 2 6zi(t)- cos 2 9zf (t) 
e i2$ gin 2 q z * ^ _ cog2 0^ ^ i — sin cos (e"^ (tj+e^zf (t) 

(26) 

Recall that in our example zi(t) = z*(t) = \ llfcLi,^ cos (2gikt). The off- 
diagonal elements of the reduced density matrix f[2ll are proportional to 
zi(t) [or zf(t)] in any basis that has off-diagonal elements, and therefore at 
least in the general spin-| model the common idea of the relation between 
off-diagonal elements of the density matrix and quantum correlations (coher- 
ence) is valid, and the coherence function E(pi,t) = 2\zi(t)\ describes well the 
time evolution of coherence. As explained in Section |2~T1 ifiV^ooAt^oo, 
then z(t), z*(t) — ► 0, and thus the diagonalisation of the reduced density ma- 
trix (in infinite model) happens in all possible basis sets. In a finite model, 
Poincare recurrences are present. The dynamical effects (decoherence and 
recoherence) do not depend on the chosen basis set, unlike in the decoherent 
histories approach, and thus, the preferred basis problem is not realised in 
density matrix formalism. 



6.2 Poincare recurrence 

The interesting thing to notice is that our quantum system (of N particles) 
is a closed and finite quantum system. Closedness results in the idealistic 
coherence being a constant of motion, and finiteness in that the realistic 
coherence S re (t) has a finite Poincare recurrence time. There is an elementary 
procedure for an upper estimate of the recurrence time of a system that 
consists of N subsystems with M = 0.5(iV 2 — N) fixed periods Tj. From the 
construct 



where rii and di are smallest possible natural numbers, the upper limit of 
Poincare recurrence time is obtained 



M 



Ip = T uni tI]>- (28) 

i 

The procedure of evaluating Tp of a particular simulation run is rather 
easy to do numerically, but with a drawback of limited accuracy. Thus, the 
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presented Poincare recurrence time analysis gives only an estimate of the 
magnitude of the recurrence times. The best fit for the data of the Figure [3] 
is given by the function 

T P {N,n p , V ,e,D) = ir^n^e 3 - 07 ^-^. (29) 

In Ref. [8], Poincare recurrence time of a similar system is proportional to 
T P ~ A!. The difference between the estimates may be a result of limited 
numerical accuracy on applying our method. Another possibility is that there 
is unintentional "double-counting". Nevertheless, our estimate gives at least 
an upper estimate for recurrence time. 

It seems clear that the system may return to its initial position, but the 
recurrence time grows fast with respect to A. For example, for A = 100, 
D = 3, e = 1, n p = 10 30 m" 3 and 77 = 8.22 x 10 43 e 2 HzmC" 2 (electro- 
magnetic interaction with charge of electron e) the estimated recurrence time 
is Tp ~ xo 13183 s and for ~ iV! the recurrence time is Tp ~ 10 142 s that is still 
a very long time compared with the age of the universe 7\j ~ 4 x 10 17 s! 



6.3 The dependence of coherence decay on system vari- 
ables 

As shown before, the decay obeys well the coherence profile = (1 — 
c ) e -(*/*d) _|_ c The floor of fluctuation level behaves as c ~ lO^ 5 if D + 



1 > y A/200exp (e/2). To derive the dependence of coherence decay on 
relevant system variables has two separable steps remaining: to study power 
parameter C and decoherence time parameter £<j- 

The power parameter C does not seem to be a function of the particle 
number A, as Figure 0] clearly shows. The function 

f(D, e) = 1.97 (l - 0.93 exp (-OMD 1 - 3 ^- 1 - 68 )) (30) 

gives the best fit for the values of C of the fitting function £(t). As both D 
and e increase, the deviation of data points in parameter C increases. 

Quite a good approximation of the functional behaviour of the decoher- 
ence time-related parameter t& of is given by 

//Us -0.085(D-l)/e 

*=<-V D (4) (») 

x(o.29e- a79 ( D - 1 ) el/4 - ai3 ^- 3 - 4 )+0.17^((D-2) 2 ) a ^ 

We assume that the accuracy of the fit Td is limited to the close environ- 
ment of the parameter space where D,e < 10, since the form of the fit is 
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very complex and most likely it does not follow the physical behaviour of 
the system accurately enough; it is only a fit on the observed behaviour. 
Moreover, the deviation in data points increases as both D and e increase. 
The most evident reason for this is the fact that N = 200 particles is still a 
small number in a short-range interaction parameter space (e > 3) when the 
number of particles pro dimension is very small (N/D < 50). In the future, 
this inconvenience should be easily overcome as the available computational 
capacity increases. The behaviour of Td with respect to the particle number 
N is more reliable because of the few data points with N = 400 and N = 800, 
and because particle number-related functions factorise out of the complex 
form in equation (l32l) . 

Despite the fact that for each parameter point there are only 100 sim- 
ulation runs and there are accuracy problems as both D and e increase, 
the behaviour of f(D, e) and Td(iV, n p , r], e, D) is quite reliable for physically 
relevant cases. Physically most interesting results are: 

• The decoherence time Td behaves as ~ 

N -const.{D-l)/e respect to 

the particle number N. Is there a constant level that decoherence time 
approaches as the particle number increases? To answer this question, 
one seems to need simulations with N that is at least one or two orders 
of magnitude greater. 

Moreover, the fit has (D — 1). Is this "—1" a physical property of the 
model or a result of numerics? We suspect the latter, but there was no 
suitable fit of Td without " — 1". 

• The upper limit of the exponent f(D,e) of the decay profile S(t) = 
exp[-(t/r d ) /(D ' e) ] is ~ 2. 

• The decoherence time Td seems to have a lower limit with respect to 
D and e. For iV = 200, it is around ~ 0.07 in units of simulation time 

r] (n p [m\Y /D t. 

• The behaviour of Td in the parameter space (D, e) is dominated by two 
different trends. One trend is a Gaussian peaking near (D = l,e = 3.4) 
and the other one is a runaway solution when both D and e are large. 
We suspect that at least a part of that Gaussian is a physical property 
of the model, but also that most of the runaway solution is caused by 
too few particles in the numerics. 

• For D = 3, N = 100 and n p = 10 30 m~ 3 , the electro-magnetic in- 
teraction (e = 1) results in decoherence time r d = 4.3 x 10~ 17 s 
and exponent f(D,e) = 1.87. For spin-spin interaction (e = 3) and 
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rj = 3.27 x 10~ 26 m 3 Hz (for 6 Li) with the same setup, the decoher- 
ence time is r d = 3.3 x 1CT 5 s with exponent f(D,e) = 0.80. For 
pure spin-spin interaction, the experimental observation of coherence 
decay seems possible in few-particle systems, but only iff particles are 
isolated enough from other (stronger) interactions, most notably from 
electro-magnetic interaction. 

While a more sophisticated study needs more computational power, ob- 
vious results of this research are that our definition of coherence is con- 
firmed by the case-study, and that the coherence decay profile is S re (i) = 
(1 — c) exp —(t/TdY + c m the Heisenberg model with r]/r £ -dependent in- 
teractions in various dimensions. 



7 Discussion 

In the past, there have been arguments that while decoherence in closed sys- 
tems is practically a working concept, closed systems do not, in principle, 
experience decoherence. A famous argument against decoherence in closed 
systems is from Bell [23], [60]. Bell's argument is mainly that if the universe 
starts in a pure state, it will always remain in a pure state, no matter how 
quickly the off-diagonal elements of the reduced density matrix decrease and 
how small they will become. Bell claims that this gives, in principle, a pos- 
sibility to make such a measurement that will show quantum interference. 
However, there is not such a measurement device, even in principle, that can 
perform the measurement on the universe. We demonstrated this with the 
Clever Chinchilla, who is a supernatural-like observer outside the universe 
(closed system). She can guess right the wave function of the universe, along 
with the interaction Hamiltonian, and then calculates the time evolution of 
the universe. In quantum measurements, one often mixes the view of the 
Clever Chinchilla and the Real Observer who studies readings of a measure- 
ment apparatus inside the universe. The Clever Chinchilla is not correlated 
with the universe, while our Real Observer is totally correlated with it, and 
thus their views about a particular quantum event may differ. The Real 
Observer can never achieve the Clever Chinchilla's knowledge about the uni- 
verse. He has no means to solve by measurements what is the whole wave 
function of the universe. He only observes projections (possibly very compli- 
cated ones) of the wave function of the universe. For a detailed discussion 
about Bell's argument, one is encouraged to study references [23l 160] . 

However, the universe is a closed system and decoherence is an observed 
phenomenon in our universe. Thus, (de)coherence should be uniquely quan- 
tifiable. In this paper, we present a unique way to quantify coherence. The 
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growth of entropy has the same dynamical origin as decoherence. In con- 
clusion, they both are equally correct descriptions, also in closed systems. 
Moreover, they grasp different aspects of a quantum system. Coherence 
as defined in this article expands knowledge of the properties of quantum 
systems beyond what is possible to obtain from entropy. First of all, coarse- 
grained coherence S re could be used to calculate the decoherence time of the 
system. It is impossible to obtain decoherence time from entropy. Second, 
it is the quantity that restricts the observation of quantum properties of a 
system more dramatically than entropy. They both give different kinds of 
constraints that a closed quantum system must obey. It is possible to have 
two systems with equal coherence but different entropy, and it is also possible 
to have two systems with equal entropy but different coherence. The latter 
group of systems is very interesting for quantum mechanical engineering - 
especially for studies related to quantum computation. Third, coherence 
may be useful in quantifying entanglement. Quantum computing still lacks 
a unique method to quantify "the usefulness of entanglement" beyond two 
entangled systems. Our proposal for mutual entanglement E = — S rc 
may be worth studying in this respect. 

The view, as the Real Observer inside the universe perceives it, can be 
obtained by using realistic coherence and entropy, which are effective theory 
descriptions about the knowledge of the Clever Chinchilla. Thus, the use of 
reduced density matrices is justified. Decoherence seems to work also in prin- 
ciple in closed quantum systems, and consistency problems can be avoided 
by studying decoherence of reduced density matrices instead of decoherent 
histories. 

To demonstrate the usefulness of coherence as defined in this article, we 
solved the time evolution of (realistic) coherence in the closed and finite spin 
model and noticed that the decay profile of coherence does not necessarily 
follow the commonly expected exponential decay. The trivial result that the 
diagonalisation of reduced density matrices occurs in all possible basis sets 
suggests that in "real" quantum experiments and measurements, observed 
noise or "measurement errors" in the behaviour of highly quantum correlated 
systems may be results of a weak, unmodelled interaction that causes deco- 
herence. Increased understanding on this topic may lead to better modelled 
quantum measurements, better measurement technology and technical solu- 
tions which may help to keep coherent systems as coherent as possible in the 
real universe. 
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Figures 




Figure 1: The behaviour of coherence function E(t) in various scenarios (a) 
compared with single particle correlations ze(t) (of randomly chosen 6 th par- 
ticle) and long-time fluctuations with system parameters N — 20, D — 1, 
e = 1, n p = 1, (b) compared with fitting function with system parame- 
ters N = 100, D — 1, e — 1, n p — 1, and (c) compared with fitting function 
with system parameters N = 100, D — 3, e=l, n p — 1. The unit of 
simulation time is T = rj (2n p [m\Y^ D t, where t is the real time. The fitting 
parameters (b) t d = 0.1212, C = 1.0070, c = 2.3 x 10~ 16 , x 2 = 0.8420 and 
(c) t d = 0.1114, C = 1.8443, c = 1.5 x 10 -21 , x 2 = 0.0744 demonstrate 
the accuracy of proposed coherence decay profile, and that fluctuations from 
theory are insig nificant (c) when D + 1 > J N/200 exp (e/2). 
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Figure 2: Expectation value of fluctuations (a) and average of weights (b) in 
respect to dimension D and the potential term e with simulation parameter 
values N = 200 and n p — l. These plots demonstrate that the proposed co- 
herence decay profile is very accurate especially if D + l > JN/200exp (e/2). 
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Figure 3: An upper estimate of Poincare recurrence times for simulation 
data (gray dots) with parameters D = 1, e = 1, and n p — 1. The unit of 
simulation time is T = rj {n p [m\) e ^ D t, where t is the real time. Best fit (solid 
line) differs considerably from ~ iV! type of behaviour (dashed line). 
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Figure 4: The dependence of power parameter C on system variables N, e 
and D. The density is n p — 1. The slices of parameter space are D = 3, 
e = 1 (a), and N = 200 (b). For certain parameter values (N, D, e) there 
exist 100 data points. It seems that there is no relation between C and the 
number of particles N (a). With small values of N, the standard deviation 
of weights is greater than average of weights, which explains deviations of 
weighted average from a constant value. With respect to dimension D and 
potential parameter e, the best fit for power parameter C is given by the 
function f(D,e) = 1.97 (1 - 0.93 exp(-0.65D L35 e - L68 )) (b). Note that for 
large D and e, the deviation in data points gets wider. 
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Figure 5: The dependence of decoherence time parameter id on system vari- 
ables N, e and D. The slices of parameter space are D = 3, e = 1 (a), D = 1, 
e = 2 (b), and N = 200 (c). The density is n p = 1 and the unit of simulation 
time is [id] = r\ (n p [m]) e ' D i, where i is the real time. For certain parameter 
values (N, D, e) there exist 100 data points. A rather good approximation 
of functional behaviour of id is given by of equation ( 1321) . Note that for 
small iV (a,b), large D and e (c), the deviation in data points gets wider. 
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